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Thin paper treats the problem of analyzing a first-come first-served 
queuing system, in equilibrium, when subjected to a peaked input {e.g., 
traffic overflowing a trunk group with Poisson input) . The basic GI/M/N 
{renewal input to N exponential servers) queuing result is used, together 
with each of two models for representing peaked traffic, the Equivalent 
Random {E-R) model and the Interrupted Poisson Process {IPP) model. 
The equilibrium virtual delay distribution is derived and compared with 
the equilibrium distribution of delays seen by arriving calls. Numerical 
examples are presented, along with comparisons of results using both the 
above models. The results show that delays can be quite sensitive to 
peakedness. 

I. INTRODUCTION 

It is well known, from analysis of blocking in trunk groups, that the 
blocking seen by peaked traffic (e.g., traffic overflowing a first-choice 
trunk group with Poisson input) can be significantly larger than block- 
ing seen by Poisson traffic with the same intensity. In this paper, we 
are interested in determining the effect peaked traffic has on delays in 
queuing systems. The analysis was motivated by a study of sender 
attachment delay in Crossbar Tandem switching machines receiving 
alternate routed (peaked) traffic. 

We treat the problem of analyzing a first-come first-served queuing 
system with peaked input for the situation where there is no idle 
server if there is a waiting customer. The basic tool is the GI/M/N 
queuing result which requires a characterization of the input process 
in terms of the Laplace-Stieltjes transform of the interarrival time 
distribution. This characterization is provided using Wilkinson's 
Equivalent Random 1 (E-R) model where the peaked input is modeled 
as an overflow process from a finite trunk group with Poisson input. 
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The size of the finite trunk group and intensity of the Poisson input are 
chosen so that the overflow process produces the desired mean and 
peakedness (variance to mean ratio of trunks up on an infinite trunk 
group). The E-R method has been widely used in analyzing blocking 
in trunking networks and is considered our basic model when the source 
of the traffic peakedness is from alternate routing. 

A second model, which gives a much simpler characterization of the 
peaked input (both computationally and analytically), is the Inter- 
rupted Poisson Process. 2 Here the peaked traffic is considered to be the 
output process of a switch, with Poisson input, where the switch is 
opened and closed for independent, exponentially distributed time 
intervals. The parameters of the switch are chosen to match either 
the first two or three moments of trunks up on an infinite trunk group 
with the corresponding moments obtained for the E-R model. 

Comparisons between results using each of the above models are 
presented along with a set of numerical results which show that delays 
can be quite sensitive to peakedness. 

For Poisson traffic, the virtual delay distribution (time congestion) * 
is identical with the delay distribution seen by arriving calls (call 
congestion). This is not the case for peaked traffic. Since in some 
applications measurements form estimates of time congestion (e.g., 
SADR measurements), it is of interest to relate the time and call 
congestion quantities. Numerical results show the sensitivity of the 
relationship to peakedness. 

ii. gi/m/n queuing results 1 

Consider a recurrent input 1 to an N server queuing system. The 
service times are independent, exponentially distributed with the mean 
service time given by 1//X2- The customers are served in their order of 
arrival, and there is no idle server if there is a waiting customer. 

Let F(x) denote the distribution function of the interarrival times. 
The mean interarrival time, 1/X 2 , is given by 



W 

A 2 Jo 



xdF(x) (1) 



* Time and call congestion are commonly used in trunking analysis (BCC system ; 
i.e., blocked calls cleared). In the delay case (BCD system), we use call congestion for 
the delays seen by arriving calls and time congestion for the virtual delay. See appendix 
for precise definitions. 

T The reader is referred to chapter 2 of Ref. 3 for a more detailed mathematical 
description of these results. 

* The peaked input will be modeled bv recurrent processes. For these processes, 

F(0+)=0. 
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and the Laplace-Stieltjes transform of F(x) is given by 
*(s) = / e-"dF(x). 

Jo 

We define the load/server (sometimes called occupancy) as 



p = 



JVV 
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(2) 



(3) 



The following result is available to us: 3 If p < 1, then the equilibrium 
delay distribution (as seen by arriving calls, i.e., call congestion) exists 
and is given by 



Prfjdelay > T~\ = 



1 - 



exp[-A r /x 2 (l - w)T]. 



(4) 



The exponential delay distribution is seen to be a function of two 
parameters, co and A. The parameter co is the solution, in (0, 1), of the 
equation 

« = *rjV M2 (l - w )]. (5)* 

For p < 1, this equation is known to have a unique solution in (0, 1); 
furthermore, the solution can be found by successively iterating on 
(5); i.e., 

o),+i = *[JV// 2 (1 — co t )] 



with coo G [0, 1). If we define 
and 



*; = *(ift) 



then the parameter A is given by 
A = = 



(6) 

(7a) 
(7b) 

(8) 



.1^^ + y?i C>(1 - * y )[JV(l -co) - j'l 

Thus, given the characterization of the input in terms of $(s), eqs. 
(6), (7), and (8) provide the means to compute the equilibrium delay 
distribution as seen by arriving calls (4) . The mean of the equilibrium 

* In Poisson traffic, the solution of (5) is w = p. In some cases, to be treated later, 
(5) can be solved in closed form. 
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delay distribution is given by 

W*l m nJ- ,)■ • (9) 

In the next two sections we discuss the Equivalent Random and 
Interrupted Poisson models for generating peaked traffic and charac- 
terizing 3>(s). 

III. EQUIVALENT RANDOM MODEL 

The E-R model treats peaked traffic as an overflow process from a 
finite trunk group with Poisson input. The holding time for the trunk 
group is assumed to be exponential with mean l//*i-* The number of 
trunks and the intensity of the Poisson traffic are chosen so the mean 
and variance to mean ratio (peakedness) of trunks up on an infinite 
overflow group closely match the desired mean and peakedness. If we 
denote the desired mean and peakedness by m d and z d respectively, 
then Rapp's formulas 4 for the Poisson load (a eq ) and number of 
trunks (c), 

a eq = ^ = m d z d + Sz d (z d - 1) (10) 

Mi 

c = aj w ' + " )-„,,-!, (11) 

q \ m d + z d - 1 / 

yield an overflow process which approximates the desired process. It 
should be noted that use of (10) and (11), or truncation of c from (11), 
will often lead to overflow traffic with mean and peakedness different 
from (but usually close to) the desired values. To quantify this effect, 
the actual mean and peakedness should be computed using 

m = a eq B(c, a eq ) (12) 

. - [i — + e + i^-J . ™ 

where B is the Erlang B function. 

Takacs 3 , Chapter 4, shows that the Laplace-Stieltjes transform of the 
interarrival time distribution of the overflow traffic is given by 

*w = t (;) ±% (. + w/je C 1 i% (« + w. w 

j=0\j/ Kq*=0 I 3-0 \ J / Aeq»=0 

* Note that we can have m ^ hi. For example, if overflow traffic from trunks is 
offered to a group of senders, l/m is of the order of minutes, whereas 1/in is of the 
order of seconds. 

tr Tne value c is considered an integer. In practice, we round c given by (11) up 
and down and choose the one that gives an actual z closest to z d . The actual m and z 
are computed from (12) and (13) using the rounded values of c. 
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where the empty product is unity. Note that <t>(s), given by (14), has 
to be repeatedly evaluated in (6) for the solution of (5) and in (8). In its 
present form, (14) is unsuitable for computation (for large c) because 
of numerical problems. In order to avoid these problems, we use a 
recursive method for the evaluation of $(s) developed by A. Descloux. 5 
If we use the notation $ c (s) to denote the dependence of <£>(s) on c, 
then $ c (s) satisfies 

[**(«) y-i = JL + 1 + ti k -£, fc **-i ( a) (15) 

Aeq Aeq A e q 



with initial condition 



*•(.) = ,-*-. (16) 



The solution is thus obtained as follows: Using (15) to evaluate 
<p(s), we iterate on (6) to find the w parameter and subsequently the A 
parameter (8). Having evaluated the w and A parameters we can 
compute Pr[delay > T~\ from (4). 

We now turn our attention to the interrupted Poisson process model 
for generating peaked traffic, the resulting simplifications, and com- 
parisons with the E-R model. 

IV. INTERRUPTED POISSON PROCESS MODEL 

This model, suggested by W. S. Hayward and analyzed by A. 
Kuczura, 2 treats peaked traffic as a Poisson process modulated by a 
random switch where the switch is opened and closed for independent, 
exponentially distributed time intervals. The importance of this process 
is that it can provide a simple and accurate approximation to overflow 
traffic. 

This model contains three parameters, the intensity of the Poisson 
process into the switch, X, calls per second; the mean open time of the 
switch, 1/w, seconds; and the mean closed time of the switch, 1/f, 
seconds. If we choose 

h = Ai = mz + 3z ( z _ 1 ) j (17) 

~ = «. = j- § Lm + 3* - 1], (18) 

and 

£-* -[£-i]«. (19) 

then the mean and variance to mean ratio of trunks up on an infinite 
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trunk group, with mean holding time 1/jui, will be m and z respectively. 2 
This corresponds to the two-parameter match of Ref. 2. The so-called 
three-parameter match is obtained by matching the first three moments 
of trunks up on an infinite trunk group with the corresponding moments 
that would be obtained using the E-R model. In this case, c and a eq are 
computed from (10) and (11) and the switch parameters are obtained 
from 



— — A a — fleq 

Ml 



82(81 — So) — 80(82 — 81) 



and 



(Si - So) - (82 - Si) 

a, _ _ 8± T A, — fl e q8l 1 
Hi A, \_ Si — 8u J' 

7. co, I A a — a eq 8o 1 

Ml °eq L do J 



where the 8*, defined in Ref. 2, are given by 

8 = B(c, a eq ) 
and 

, C ~r K fleq 1 Aeq s 



k 



k 



(20) 
(21) 

(22) 

(23)' 
(24)t 



For a given mean m and peakedness 2, it has been shown 2 that the 
Laplace-Stieltjes transform of the interarrival time distribution of the 
output process of the switch is given by 



kiTim k 2 r 2 fn 

$(s) = — r— 



s + rim ' 8 + r 2 fj.i 
where the parameters r h r 2) k 1} and k 2 are given by 



(25) 



ri = i[A. + a>, + 7. + V(A, + w. + t.) 2 - 4ii.w.] (26) 
r 2 = HA. + a,, + 7. - V(A, + co. + 7*) 2 - 4A a co„] (27) 
A, — r 2 



ri - r 2 
fc 2 = 1 — ki 



(28) 

(29) 

and A s , co., and 7. are given by (17), (18), and (19) for the two-param- 
eter match and by (20), (21), and (22) for the three-parameter match. 

* B is the Erlang B function. 

t This equation can be simply obtained from equation (1.15) in the appendix to 
Ref. 1. 
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With (25) defining $(s), eq. (5) becomes a cubic in the w parameter 
of the delay distribution. Dividing the cubic by the known root at unity 
gives the desired root in (0, 1) 



/ 1+flrN. \ 
"=( 2 )- 



with 
and 



Ul +at + a»y 


NltoV + AW 


rim 




T2fil 

a2 = W>' 





(30) 



(31) 



(32) 



Thus a:, given by (30), together with (8) specify the equilibrium delay 
distribution. Note that the iteration procedure (6) has been eliminated 
and that 4> ; , defined by (7a), is simple to compute using (25). 



V. TIME CONGESTION 

For Poisson traffic, the virtual delay distribution (time congestion) 
is identical with the delay distribution as seen by arriving calls (call 
congestion). This is not the case for peaked traffic. Since in some 
applications measurements form estimates of time congestion (e.g., 
SADR measurements), it is of interest to relate the time and call 
congestion quantities. 

We define time congestion as the delay in being serviced experienced 
by a fictitious call arriving at an arbitrary time t when the system is in 
equilibrium. * It is shown in the appendix that, in equilibrium, the rela- 
tionship between call congestion (CC) and time congestion (TC) is 
given by 

Pr[TC > Tl = ? PrlCC > 7 1 ] (33)* 

if p < 1 and the interarrival time distribution, F(x), is not a lattice 
distribution.* Pr[CC > 7 1 ] is the delay distribution seen by arriving 
calls and is given by (4). 

VI. NUMERICAL RESULTS AND DISCUSSION 

One measure of system performance of possible interest is the mean 
delay experienced by arriving calls, given by (9), versus CCS offered 

* See appendix for a more precise definition of time congestion. 

* For Poisson traffic, w = p giving Pr[TC > T] = Pr[CC > T]. 

1 The E-R model and Interrupted Poisson model clearly satisfy the nonlattice 
hypothesis. 
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to the servers where 



#CCS = SG^m 



with m given by (12) in the E-R model and 

m = A, 



+ 7. 



(34) 



(35) 



for the interrupted Poisson model. Here A a , y„ and co. are given by 
(17), (18), and (19) for the two-parameter match and by (20), (21), 
and (22) for the three-parameter match. Figures 1(a) and 1(b) are plots 
of the mean delay characteristics for each of the three models of 
interest. The values of peakedness z ranging from 1 to 3 are presented 
on Fig. 1(a), and the z = 4 results are plotted on Fig. 1(b). The param- 
eters of this example are N = 18, l//*i = 180 seconds, and 1/m = 7.6 
seconds. 

It is seen that, while the two-parameter results tend to overestimate 
both the E-R and three-parameter results, the differences are indis- 
tinguishable (for the entire CCS range shown) up to 2 = 1.5 and small 
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Fig. 1 — Mean delay characteristics. 
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Fig. 2 — Delay distribution. 

up to about 2 = 2. The three-parameter results* are seen to be in- 
distinguishable from the E-R results up to z = 2.0 and close up to 
about z = 3. In all cases, the E-R results tend to lie between the two- 
and three-parameter results. We see in Fig. 1 (b) that the results differ 
greatly for z = 4. The results of Ref . 6 show that fixing the equivalent- 
random mean and variance for a renewal process does not necessarily 
tightly tie down the blocking in the BCC case. We are observing the 
same phenomenon here. In fact, nonnegligible discrepancies between 
the E-R and three-parameter results for larger z's can occur despite 
a matching of three moments. 

Specific delay distributions are shown in Figs. 2(a) and 2(b). An 
observed property of the results are that the o> parameter (5) for the 
two-parameter case exceeded the oj parameter of the E-R model which, 
in turn, exceeded the w parameter for the three-parameter case. This 
explains the slope differences. This also tends to order the T = 
results as shown. In many cases, the A parameter of the distribution 

* In order to avoid severe numerical problems in computing the three-switch 
parameters, double precision was used in eqs. (20) through (24). 
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(8) for the two-moment match exceeded the A parameter for the E-R 
model which, in turn, exceeded the A parameter for the three-param- 
eter match. From these figures, we again observe the closeness of 
results for z = 2 and the large differences for 2 = 4. 

Figures 3(a) and 3(b) plot the load service relationship P(delay > 2.5 
seconds) versus offered load for each of the three models. The compari- 
sons again exhibit the same characteristics that were seen for the mean 
delay results [Figs. 1(a) and 1(b)]. 

These results show the extreme sensitivity of the queuing system 
performance to the peakedness of the input process. They also give some 
insight into the region where two- and three-parameter results are 
expected to be closest, i.e., low peakedness and high congestion. When 
seeing the discrepancies between the two- and three-parameter matches 
and the E-R model, we may question which is the bench mark. If the 
peakedness arises from alternate routing, the E-R model seems basic 
since it is an overflow model which has been shown to accurately 
describe superposition of overflows. 1 This has led to its wide use in 
analyzing blocking in trunking networks. It should also be noted that 
the parameters of the interrupted Poisson process have been chosen to 
match the moments of the E-R model. 
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Fig. 3 — Load service characteristics. 
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Fig. 4 — Call and time congestion. 



Although there are discrepancies between the two- and three- 
moment matches and the E-R model for high z's, the IPP is a close 
approximation to the E-R model for a wide range of practical z's. 
Furthermore, it should be emphasized that it provides a convenient 
method of analysis using birth and death equations (and simplified 
Laplace-Stieltjes transform) in many cases where the E-R model is 
intractable. 

Figures 4 and 5 result from applying (33) to the example under 
consideration using the E-R model. Figure 4 shows the load service 
relationship for both call and time congestion. The call congestion 
results, taken from Figs. 3(a) and 3(b), are reproduced here for com- 
parison. It is seen that the time congestion (TC) results consistently 
fall below the call congestion (CC) results.* Figure 5 shows the CC- 
to-TC ratio (w/p) as a function of peakedness. While for a given load 
the time and call congestion can differ substantially, the decrease in 
load that makes up the difference may be relatively small. 



* It has been shown by R. P. Marzec that for the Interrupted Poisson Pro- 
cess, Pr[TC > 0] < PrlCC > 0]. This implies p/co ^ 1, which in turn implies 
Pr[TC > T] £ PrlCC > T]. 
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Fig. 5 — Call-to-time-congestion ratio. 
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APPENDIX 

Time Congestion 

Let £(0 denote the state of the system (number of customers waiting 
or being served at time t) and let P k (t) = Pr[£(t) = fc]. Moreover, 
define £„ = |(C) (i.e., the number of customers in the system just 
prior to the arrival of the nth customer). Takacs shows (Ref. 3, 
Theorem 1, p. 148) if p < 1 then the limiting distribution 

Um P[£„ = *] = P k (ft-0,1, -••) (36) 

n-»oo 

exists and is independent of the initial distribution. Furthermore, he 
shows that 

P k = Au> k ~ N (k « N, N + 1, • • •), (37) 

where A is given by (8) and a> by (5). It is also true that the above 
holds for k = N — 1, i.e., 

P k = Ao> k ~ N (k = N - 1, N, N + 1, • • ■)• (38) 

If we denote by 77 „ the waiting time of the nth customer, then the 
equilibrium delay distribution, 

W(x) = lim W n (x) = lim P[> n ^ a;], 
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exists and is given by 

w(x) - % p ' + ,4 p 'f '"" ?^n ***• (39) 

which reduces to (4). We have 

w(x) = L 1 Pi + t PMx), 

where Ij(x) represents the integral in (39). The complementary dis- 
tribution, W(x), is given by 

W(x) = 1 - W(x) = £ P£l ~ Ij(x)l (40) 

If p < 1 and the interarrival distribution F(x) is nonlattice, then the 
equilibrium time congestion probabilities given by 

P t * = limP*(0 (k = 0, 1, ••■) (41) 

exist and are independent of the initial state. Furthermore, 

P* t = pP*_, (fc = N, N + 1, • • •)■ (42)* 

At this point, we are able to evaluate the equilibrium time congestion 
distribution under the hypothesis p < 1, and F{x) is not a lattice 
distribution. 

Denote y\ (t) as the waiting time of a fictitious arrival at time t and let 

W(t,x) = PrZ v (t) =£*]. (43) 

We have 

W(t, x) = N ± l Pfl) + ± Pj(t) I' e-^ ^^ N^dy (44) 

>=0 j=N JO {J — I\) \ 

since v (t) = if £(0 < N. And if £(0 = j ^ N, then fictitious arrival 
must wait for j + 1 — N successive departures. These departures 
follow a Poisson process with intensity JV/x 2 . Taking the limit and using 
(41) we have 

W*(x) = Urn W(t, x) = N £p* + £ P)lj(x). (45) 

Using (42), the complementary distribution is given by 

W*(X) = P t iVi[l - //(*)]. (46) 



* See Theorem 2, p. 153 of Ref. 3. 



1228 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1973 

From (38), (40), and (46): 

W*(x) = ?W(x)- (47)* 

CO 

This result is given by Ref. 7, p. 229, for a single-server queue; the 
multiserver case is given here for completeness. 
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* Note that, for Poisson traffic, w = p, which gives W*(x) = W(x). 



